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ABSTRACT 

Recent EUV spectroscopic observations indicate that slow magnetosonic 
waves are present in active region (AR) loops. Some of the spectral data were 
also interpreted as evidence of fast (~ 100 — 300 km s _1 ) quasi-periodic flows. We 
have performed three-dimensional magnetohydrodynamic (3D MHD) modeling 
of a bipolar AR that contains impulsively generated waves and flows in coronal 
loops. The model AR is initiated with a dipole magnetic field and gravitationally 
stratified density, with an upflow driven steadily or periodically in localized re- 
gions at the footpoints of magnetic loops. The resulting flows along the magnetic 
field lines of the AR produce higher density loops compared to the surround- 
ing plasma by injection of material into the flux-tubes and the establishment of 
siphon flow. We find that the impulsive onset of flows with subsonic speeds result 
in the excitation of damped slow magnetosonic waves that propagate along the 
loops and coupled nonlinearly driven fast mode waves. The phase speed of the 
slow magnetosonic waves is close to the coronal sound speed. When the ampli- 
tude of the driving pulses is increased we find that slow shock-like wave trains 
are produced. When the upflows are driven periodically, undamped oscillations 
are produced with periods determined by the periodicity of the upflows. Based 
on the results of the 3D MHD model we suggest that the observed slow mag- 
netosonic waves and persistent upflows may be produced by the same impulsive 
events at the bases of ARs. 

Subject headings: Sun: activity - Sun: corona - Sun: flares - Sun: oscillations - Sun: 
UV radiation: waves: magnetohydrodynamics (MHD) 
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Introduction 



Observations of slow mode oscil 
in X-ray emission in the past (e.g. 



ations in hot loops associated with flows were seen 



Svestka 



19941 ) . Observations with the Solar and 



Heliospheric Observatory (S OHO) Solar Ultra v iolet Measuremen t s of Em itted Radiation 



(SUMER) instrument (e.g., 



Kliem et al. 



2002 



Wang et al. 



2002 



2003b 



aj) show that 



standing slow magnetosonic waves are excited in hot (~ 6 — 9MK) coronal loops by 
impulsive flare-like phenomena. Recent Hinode satellite Solar Optical Telescope (SOT) 
observations, and Transition Region and Coronal Explorer (TRACE) satellite Extreme 
Ultraviolet (EUV) imaging observations show that some coron al loop oscillations are 



associated with 



Ofman fc Wang 



lows and jets in cool (~ 1MK) coronal loops ( lAschwanden et al 



2002 



20081 ). The difficulties in distinguishing between slow magnetosonic waves 



and quasi-periodic flows in some EUV time-sequence images of coronal loops are well known 



(e.g. 



Berghmans fc Clette 



19991 ) . The use of spectral data in addition to stereoscopic 



EUV imaging to help resolve whether the observed event is slow magnetosonic wave 



or qu asi-periodic flow is a promising approach (e.g., 



Marsh et al 



2009 



Marsh fe Walsh 



20091 ). However, the quasi-periodic signatures propagating along the loops seen in EUV 
difference images and quasi-periodic Doppler blue shift in spectral emission lines are still 
being debated regarding their wave/flow nature due to the limited quality of the present 



(S/N) ratios ( 


De Pontieu & Mcintosh 


2010: 


Verwichte et al. 


2010: 


Wane: et al. 


2009. 


2010 


Tian et al. 


2011 


)■ 



Linear MHD models of slow magnetosonic waves in coronal loops with steady 
background flows devel oped in the past were limite d to straight cylindrical structures, 



and linear waves (e.g., 



Nakariakov fc Roberts 



19951 ) . Impulsive excitation of slow 



magnetosonic waves without flows in bipolar magnetic structures was recently considered 
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in 2D MHD ( 


Sc 


wa et al. 


2005 


Selwa & Ofman 


2009 


) models. 



wa et alJl2005L 120071 ). and recently in 3D MHD f lOfman fe Selwal 12009 



20091 ) models. In these studies it was found that the simultaneously excited 



fast magnetosonic waves could facilitate the establis hment of standing s l ow m agnetosonic 



Schmidt fe Ofmanl (120111 ) studied 



waves in curved magnetic field geometry. Recently, 
slow magnetosonic waves in an expanding coronal loop associated with a Coronal Mass 
Ejection (CME) using 3D MHD. The propagation and dissipations of slow magnetosonic 



waves in coronal p 



( lOfman et al. 



umes with background solar wind flow were studied using 2.5D MHD 



20001 ) . Oscillations of multiple cylindrical coronal loop strands that include 



coupled slow waves, transverse waves and the nonlinea r interaction with background steady 



flow were studied using 3D MHD model (lOfman 



20091 ). 



The present study expands previous models by including for the first time impulsive 
and non-steady flow injection that excites the waves in magnetic loops of a bi-polar active 
region (AR) in the full 3D MHD model. The nonlinear interactions of the oscillations and 
flows are studied in realistic AR geometry in hot coronal AR plasma model stratified by 
solar gravity. The results demonstrate that injection of steady and non-steady flows at 
the footpoints of coronal loops excite compressional waves in longitudinal and transverse 
directions to the loop's magnetic field. The model shows that the same impulsive process 
at the corona-transition region boundary can produce flows and waves in coronal AR loops. 



2. Observational Motivation 



It was suggested that the standing s 



small (or micro-) flares at one footpoint (jWang et al 



ow-mode waves in 



rot coronal 



2003a . 



2005 



Wang 



oops 



201ll ). 



ar e excited by 



Wang et al. 



(120051 ) examined the evolution of Fe xix and Fe xxi line profiles in the initial phase for 54 
oscillations, and found that nearly half show the presence of two spectral components: the 
stationary component, and the shifted component that reaches maximal Doppler velocity 
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on the order of 100—300 km s 1 . These features suggest that the initiation of waves is 
closely associated with hot plasma upflows produced in flares. 

The oscillation event, occurring at 11:20 U T on 2002 April 17 above the north-west 



limb, is su ch a clear example (see Figures 1-4 in 



st udied b y 



by 



Wang et al. 



Wang et al. 



( 120051 )). This event was first 



( 120051 ) using SOHO/SUMER observations, and was further analyzed 



Wangi (1201 ll ) using Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) 



data to explore the trigger mechanism. The SUMER observed two brightenings in the Fe 
xix line located at the intersections of the slit with a coronal loop. The evolution of Fe 
xix spectral profiles shows the emergence of a second component, which has the Doppler 
velocity up to 230 km s -1 , measured by double Gaussian fits. The slow-mode waves were 
seen excited following this initial hot flow injection. The RHESSI observed a hard X-ray 
source near the northern footpoint of the oscillating loop, suggesting that a hot flow pulse 
produced by a microflare (e.g., via the explosive chromospheric evaporation) moves upward 
along the loop, leading to the two brightenings as observed by SUMER. This example 
indicates that the slow-mode wave may be driven by the hot flow that is produced by 
energy release near the loop's footpoint. 

Similarly, the association of continuous upflows and propagating slow waves observed in 
coronal loops near the edge of ARs by Hinode/EIS suggests that they may res ult from the 



same mechanism, associated w ith tiny nanoflare-like small-scale energy release (IWang et al. 



2010 



Nishizuka fe Hara 



201 ll ). The SUMER oscillation events have high recurrence rate 



(2-3 times within two hours), consistent with the hypothesis of their microflaring origin. 
We speculate that the nanoflares could produce nearly continuous upflows near the loop's 
footpoints because of their much higher recurrence rate, while larger, lower-frequency flares 
are responsible for the discernible quasi-periodic wave disturbances that propagate upwards 
more than 70 Mm above the footpoints. A noticeable difference is detected between the 
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upflow related phenomena observed by EIS and SUMER. The EIS continuous upflows and 
propagating slow waves are most clearly observed in warm (1-2 MK) coronal loops, while 
those detected by SUMER are observed in hot (>6 MK) flaring loops. Th is observational 



difference can be understood by recent nanoflare models (see the review by iKlimchukl 120091 ) 
that predict the presence of small amount of very hot (> 5 MK) plasma. There are two 
main reasons that make hot spectral lines hard to detect: the shorter lifetime and the lower 
densities of the hot material compared the the warm plasma. Therefore, we model these 
phenomena in hot coronal loops as a first reasonable step with the aim of understanding 
the relationship between flows and waves excitations by impulsive events at the loops' 
footpoints. 



3. Numerical Model 

Here, we describe the numerical 3D MHD model used to study the hot AR loops with 
flows and waves. The resistive 3D MHD equations are solved with gravity and isothermal 
energy equation on a Cartesian 258 3 grid using the modified Lax-Wendroff method with 4th 
order stabilization (smoothing) term. The initial magnetic field is a dipole (see Figure H]) 
with gravitationally stratified normalized density given by 

P = poe [l/(10 + -^n)-0.1]/H (1) 

where H = 2kBT Q R s /(lOGM s m p ) is the normalized gravitational scale height, R s is the 
solar radius, ks is Boltzmann's constant, To the background (isothermal) temperature, G 
is the universal gravitational constant, M s is the solar mass, and m p is the proton mass, 
and po — 1. The following normalization parameters were used in the present study: 
distances in units of a = R s /10, magnetic filed B = 100 G, T = 6.3 MK (to model the 
hot coronal AR loops), number density n = 1.38 x 10 9 cm ~ 3 , resulting in Alfven speed 
Vao = -Bo/(4vrno) 1//2 = 5872 km s _1 used to normalize the velocities, the Alfven time 
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ta = cl/Vao = 11.9 s used for time normalization, the isothermal sound speed C s 



323 



km s 1 = 0.055V^4 , and the corresponding plasma f3 



2C? 



0.00605. The model is 



based on the 3D MHD model of a bipolar AR developed initially by 



Ofman fc Thompson 



( 20021 ). This model and its extensions and variants were used successfully in a number 



of studies of waves in ARs a nd in coronal loops (e.g. 



2008 



Selwa fc Ofman 



2009 



2010 



Selwa et al. 



Ofman 



2007 



McLaughlin &: Ofman 



2011al lbl). As in previous studies we set 



the Lundquist number to S = 10 4 , i.e., the resistive diffusion time scale is 4 orders of 
magnitude longer than the Alfven time. Since we are dealing with time scales O(100r^) - 
the resistive dissipation is insignificant. Numerical convergence study has shown that the 
results reported below are not affected by numerical dissipation. 

In Figure |2]the initial density, p, the fast magnetosonic speed Vf = (Vj + C 2 ) 1//2 , where 
Va = B(x, y, z)/p(x, y, z) 1 / 2 Vao is the local Alfven speed, and the plasma (3 in the xz plane 
at y = of the model AR are shown. The gravitational stratification of the density is 
evident, and the increase of Vf in the lower central region due to the bipolar magnetic field 
of the AR is apparent. The plasma (3 increases with height due to the rapid decrease of 
the dipole magnetic pressure ~ B 2 compared to the decrease of the thermal pressure with 
height. 

Nonsteady and impulsive flow injection along the field is introduced at the lower 
boundary of the model AR as 



V = V (x, y,z = z min , t)B/\B\ 



(2) 



where 



V (x, y,z = z min , t) = A v (t)V A exp 



. 2 / \ 2 

x- x \ , fy-y 



w 



with 



Wq 



(3) 



A v (t) — A , t < t < t max , 



(4) 
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for impulsively excited continuous flow or 



A v (t) = A (l - cosut)/2, t>t c 



(5) 



for periodic upflow. Here, the parameters are z min — 1, xq — 0.8 for the short magnetic loop 
(length = 2.13), x = 1.2 for the long magnetic loop (length = 4.53), y = 0, w = 0.12, 
u = 0.21, to = 0, and the values A = 0.01, and 0.05 that corresponds to subsonic upflow 
were used. The additional power of two on the RHS of Equation fl3]) models an upflow 
with sharper than Gaussian cross-sectional profile. The upflow is imposed at the lower 
boundary z = z m i n = 1 at a circular area with twice the radius 2w$ centered at the right 
footpoint of the loops. Downflows through the lower boundary are allowed at the other 



footpoint of the loops. The rest of the boundary conditions at z — 1 are: v x 







outside the upflow/downflow regions, with fixed B, density p and v z extrapolated from the 
interior points. These boundary conditions approximate the line-tied b oundary conditions 



for coronal loops (e.g., 



Lionello et al. 



1998 



McLaughlin &: Ofman 



20081 ). modified allowing 



inflow and outflow. The line-tied boundary conditions approximate the effects of the dense 
photosphere on the loop's magnetic footpoints, and ensure wave reflection necessary for 
establishment of a standing wave with nodes at the footpoints. Open boundary conditions 
are used at the top and side planes of the computational box allowing wave propagation 
outward through the boundary with negligible reflection. 

The magnetic field lines of the loops and the dependence of the loops parameters 
on the coordinate along each loop, s, are shown in Figure El It is interesting to note 
that the normalized magnetic field strength, B and the fast magnetosnic speed Vf along 
the loop decrease with height by a factor of 3 for the short loop and by an order of 
magnitude for the long loop. The significant variation of the magnetic field and the fast 



ma 



gnetosonic speed is usua 



by 



Nakariakov & Verwichte 



y neg lected in linear coronal seismology models (see the review 



20051 ) . The density along the loops varies by ~ 10% for the 
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short loop, and ~ 23% for the long loop. The small variation of the density with height 
(compared to cool 1MK loops) is the result of the high temperature (6.3MK) and the low 
height of these loops. 

4. Numerical Results 

In Figures l4lfT4l the results of the 3D MHD model of steady and periodic inflows 
and associated waves are shown for the initial state and parameters described in the 
previous section. The evolution of the flows and oscillation is also evident in the animations 
accompanying the online version of the journal. First, we report the results for the steady 
inflow cases. In Figure H] the snapshots of the density and the velocity in the xz plane at the 
end of the run (t = 280ta) of steady inflow with Vq = 0.05 are shown. The intensity scale 
indicates the magnitudes of the density and the velocity, and the arrows in the velocity 
images indicate the direction of the flow for velocities greater than 6% of the maximal value. 
The short loop results are shown in Figure H^-b, and the long loop results are in Figure Hb-d. 
It is evident that the steady inflow leads to the formation of a steady short loop in the 
low-/? region of the model AR. The short loop is filled with inflowing coronal plasma and 
reaches higher (nonuniform) density than the surrounding plasma by a factor ~ 1.5. Since 
the EUV line emission in optically thin coronal plasma is proportional to n 2 this loop would 
appear twice brighter than the same loop preceding the inflow in background-subtracted 
EUV images. In Figure [6]the cut in the xy plane at z — 1.26 of the density and the velocity 
of the AR with the short loop are shown. The locations of the upflow and downflow that 
correspond to the loop with higher than the surroundings density are evident. In addition 
the circular structures in the density and the velocity are formed due to the compression 
and rarefaction produced by the flow in the AR and the fast magnetosonic waves emitted 
in the initial impulsive stage of the loop formation. 
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The time-distance plots shown in Figure [5] for the boxed regions E-H, (Figure H]) 
demonstrate the details of the loops oscillations. The widths of the boxes are 7 gird points 
in space and the time-distance plots are smoothed with 15 points running average. The 
snapshots of the boxes are taken every 5ta- The time-distance plot E is aligned along the 
z direction centered at x = through the apex of the short loop. It is evident that the 
temporal oscillation of the loop in the z direction is small, and the density variations are 
due primarily to compressive density oscillation, following an initial transient stage (the 
time for the inflow to reach the apex of the loop). The time-distance plot F is aligned along 
the x direction centered at height z = 1.55 through the apex of the loop and shows the 
oscillatory displacement of the loop top, with evidence of compressive density oscillations. 
By comparing time-distance plots E and F it is evident that the dominant transverse (kink) 
oscillation of the loop is in the x direction (see, Figure [8] below) . The time-distance plots 
G (centered at z — 1.27) and H (centered at z — 1.82), are cuts in the x direction through 
the footpoints of the short and long loops, respectively. Note, that the positive values of x 
correspond to the right side of the plots in Figure HI Here, the transverse oscillations of 
the footpoints, as well as the effects of density compression are evident. The displacements 
of the loops' two footpoints are in-phase, i.e., the loop sways left and right in the loops' 
curvature plane, whereas the intensity oscillations near the two footpoints are in anti-phase, 
a signature of the fundamental slow mode wave, which has the density anti-nodes at the 
footpoints. These l oop oscillations are similar to the second harmonic of the vertical kink 



mode discussed by 



Wang et al 



( 120081 ) in connection with TRACE EUV coronal active 



region loop oscillations. 

Figure [7] shows the temporal evolution of the velocity components, the perturbed 
magnetic field components, and the perturbed density at the right footpoint of the short 
loop (point A = (0.59, —0.01, 1.27)) with steady inflow. At this location the velocity along 
the loop has two components in the Cartesian frame (x and z) as evident in Figure [3^. It is 
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evident that it takes about 15ta for the perturbation along the loop to reach this position. 
The loop contains both, compressive and transverse oscillations, as evident in the time 
dependence of the velocity and magnetic field component oscillations, as well as the density 
fluctuations. The transverse waves are excited by the momentum of the injected velocity 
pulse and the centrifugal force exerted by the inflow on the curved (dipole) field lines, 
with the Lorentz force acting as the restoring force. The slow mode wave period is 65ta 
(determined from the last oscillation period of An), and the oscillations damp in several 
periods due to leakage at the footpoints and outside the loop as a result of finite (3 allowing 
oblique propagation and leakage of the slow magnetosonic waves in curved magnetic 
geometry. The quarter period phase shift between the density and the velocity components 
along the loop indicate the presence of standing slow-mode waves is established quickly. 
The presence of the slow mode wave is also evident from nearly anti-phase relation between 
\AB\ = {AB 2 X + AB 2 + AB 2 z f/ 2 and An (see the vertical red dashes lines in Figure [7]) . The 
anti-phase relation is not exact due to the effects of nonlinearity, and non-modality (i.e., 
the oscillation are not exact normal modes of the loop) of the oscillations. The effects of 
nonlinearity are evident in the non-sinusoidality of the oscillations at the initial stages of 
the evolution. The calculations were repeated with smaller velocity amplitude (Vo = 0.01) 
confirming that these effects are due to nonlinearity and the magnitude of the initial pulse. 

Figure [8] shows the temporal evolution of the velocity components, the perturbed 
magnetic field components, and the perturbed density at the apex of the short loop (point 
B = (0.01, —0.01, 1.60)) with steady inflow. The red dashed line marks three consecutive 
peaks of the density perturbation, An. It is evident that AB X is in phase with An indicating 
the presence of fast mode wave and the corresponding kink oscillation of the curved loop in 
the xz plane induced at the impulsive initial stage of the inflow. The variations in An and 
AB X are caused by the oscillations of the nonuniform loop with respect to the fixed point 
B, with contribution from compressive effects (see, Figure [5]). The period of the oscillation 
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is 35ta based on the second and third peaks. This period is shorter than slow wave period 
and lon ger than the kink m ode period of the fundamental mode for a cylinder of length 



L fe.£ 



Roberts et al 



119841) given by 2L/C k , where C k = V A (2/(1 + pj ' Pi )V* = 1.1V A , 
with the ratio of the density inside and outside the loop Pi/p = 1.25/0.78 = 1.5 from the 
model calculations. Noting that Va decrease with height in the loop we get C k = 0.35Vao 
at the apex of the loop. However, the theoretical kink mode period of 12t a is a factor of 3 
shorter than the period found in the model. Thus, the observed fast mode is not a linear 
fundamental kink mode of the loop but nonlinearly driven fast mode oscillation with similar 
properties to the second harmonic of the vertical kink mode. Further justification of this 
interpretation is obtained by inspecting the phase and period relation between AB X , An, 
and AV X . It is evident that AV X period is twice as long as the AB X , An, and peaks, nodes, 
and minima of AV X all correspond to peaks of An and AB X , suggesting that the fast mode 
is driven nonlinearly by a term that depends on (AV X ) 2 . Such second order term appears in 
the x component of the momentum equation ~ P^4^f = • The nonlinearly driven 

fast mode waves damp as the driving amplitude of the slow waves decreases with time. 

Figure [9] shows the temporal evolution of the velocity components, magnetic field 
components, and the perturbed density at the right footpoint of the long loop (point 
C = (1.08, —0.01, 1.8)) with steady inflow. The long loop exhibits similar evolution and 
contains both, longitudinal and fast mode waves induced by the flow. The slow wave period 
is ~ 136ta, determined from the last half period of An oscillation. The presence of the slow 
mode wave is also evident from nearly anti-phase relation between \AB\ and An (see the 
vertical red dashes lines in Figure |H]). The anti-phase relation is not exact due to the effects 
of nonlinearity, and non- modality (i.e., the oscillation are not exact normal modes of the 
loop) of the oscillations. As in the short loop case the waves damp due to leakage of the 
slow wave outside the loop. Evidence for the formation of slow standing wave is seen in the 
nearly quarter period phase shift between the velocity along the loop and An. However, 
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in the case of the long loop the periods of the oscillations is longer as expected, and the 
formation of the slow standing wave takes longer time than in the short loop. By comparing 
the short and long loop oscillations it is evident that the periods of the oscillations and 
the wavelengths are related to the loop lengths. The ratio of their periods is about a 
factor of two, close the loop's length ratio of 2.15. Wave period values in both loops are 
in qualitative agreement with the classical expression 2 L/C* for the fundam ental mode of 



Roberts et al 



19841 ) with the short 



the slow magnetosonic wave in a straight cylinder (e.g. 
loop in better agreement. However, quantitatively the slow wave periods is smaller by 16% 
for the short loop and by 17% for the long loop case. The difference is mainly due to the 
curved magnetic geometry and nonuniformity of the loops. 

In the following, we report the results for the periodic flow injection cases. The density 
and velocity distributions in the xz plane (at y = —0.0136) cuts of the AR produced 
by periodic flow injection at the lower boundary at (0.8, 0, 1) at times t = 22.5, 105, 
188 ta for the short loop are shown in Figure [10J The loop density structure shows 
time-dependent variation due to the combined effects of the periodic flow injection and the 
compression of AR plasma by the velocity pulse. The velocity magnitude and the direction 
vary approximately in-phase with the density showing longitudinal variations along the 
expanding magnetic field of the loop. The temporal evolution of the velocity components, 
perturbed magnetic field components, and perturbed density are shown in the short loop 
at the right footpoint at point A = (0.59, —0.01, 1.27) for the periodic flow injection in 
Figure [UJ It is evident that after the short initial phase as the perturbations reach the above 
position, the loop's magnetic field, velocity, and density exhibits quasi-periodic oscillations 
induced by the periodic flow injection at the footpoint. The velocity fluctuations along 
the loop are in phase with the density perturbation, and in anti-phase with the magnetic 
field components. The temporal evolution of \AB\ shows two dominant frequencies due 
to the combined effect of the driver and the slow mode wave in the loop. No temporal 
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damping or growth of the oscillations is evident in the driven case at the steady state, since 
the fluctuation are injected at the right footpoint continuously, and exit the loop at the 
left footpoint boundary. The non-sinusoidality of the oscillations is evident and is due to 
the nonlinear interaction between the flow periodicity, the resulting variation of the local 
loop Alfven speed that affects the properties of the Alfvenic fluctuations, and the waves 
produced in the loop. The nonlinear effects are reduced when the same calculation is 
repeated with smaller injected velocity amplitude (Vo = 0.01, not shown). 

In Figure [12] the cuts in the xz plane of the AR's density and velocity injected with 
periodic flow at the lower boundary at (1.2,0, 1) at times t = 22.5, 105, 188 ta for the 
long loop are shown. It is evident that in this case the flow is divergent as expected in the 
bipolar magnetic field geometry. In the case of the long magnetic loop the density structure 
is more complex than in the case of the short loop that is entirely in the low-/? region of 
the AR as shown above. Since the long loop reaches higher altitude than the short loop, 
nonlinear steepening of the density perturbation fronts are produced due to the decrease 
of the gravitationally stratified density and the increase of the plasma /3 with height (see 
Figure [2]). Coupled slow and fast magnetosonic shock-like waves are produced as evident in 
the compression fronts in the images, and their direction of propagation that is oblique to 
the direction of the magnetic field. Here as well no temporal damping of the oscillations 
is evident. Figure [13] shows the temporal evolution near the apex of the long magnetic 
loop (point D = (0.01, —0.01,2.5)) for the periodically driven upflow case. The shock-like 
structures are evident in the temporal evolution of the velocity components, the density, 
and the magnetic field exhibiting strong nonlinear steepening. Here as well, the temporal 
evolution of \AB\ shows two dominant frequencies due to the combined effect of the driver 
and the slow mode wave in the loop. The structures are shock-like (more diffuse compared 
to shocks), since they travel less than a single scale height, and their amplitude is significant 
but not large compared to the sound speed and fast magnetosonic speed. Convergence test 
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shows that numerical dissipation does not affect significantly the shock-like waves in the 
present calculation. The nonlinear steepening of the flow fronts in the loop is due to inflow 
amplitude comparable to the sound speed at the footpoint, and the density stratification in 
long loops (see Figure [3]) . 

Finally, we summarize the results of the four cases discussed above using 3D isocontours 
snapshots of the perturbed density in Figure [TH The four cases discussed above are: (a) 
steady inflow in the short magnetic loop; (b) periodic inflow in the short magnetic loop; 
(c) steady inflow in the long magnetic loop; (d) periodic inflow in the long magnetic loop. 
By comparing between the isocontour snapshots and animations the effects of the magnetic 
loop size and the types of inflow become apparent. A steady loop is formed following 
steady inflow and induced damped oscillations in the initial transitory stage (longitudinal 
and transverse) in the short loop that lays entirely in the low-/? region of the AR. When 
the periodic flow is injected in such magnetic loop the density structure in the loop is 
dominated by the structure of the flow with small effects due to waves in the loop. In the 
long loop the steady state is not reached, and the density of the loop is variable throughout 
the evolution for both types of inflow. When the inflow is periodic the effects of the waves 
induced by the flow are significant in the long loop (for the magnitude of the flow in the 
present study) and the density structure exhibits the coupled fast magnetosonic and the 
slow wave shock-like structure in the high-/3 region of the model AR. 



5. Discussion and Conclusions 

High resolution and high cadence observations indicate that some coronal loop 
oscillations are associated with flows along the loops and with jets at loops' footpoints. 
Whether propagating intensity disturbances seen in EUV imaging observations of the 
coronal loop structures are flows or waves is still under intensive debate. On one hand 
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the propagating disturbances have been interpreted as signatures of slow magnetoacoustic 
waves since their propagating speeds are close to the sound speed in the corona, and 3 
min or 5 min periodicities suggest their association with the l eakage of the photospheric 



p- modes (e.g., 



Nakariakov et al. 



2000; 



De Moortel et al. 



20021 ). Combined stereoscopic 



EUV and spectral observations provide further evidence for slow magnetosonic waves in 



coronal loops (e.g.. iMarsh et al. 



2009 



Marsh fc Walsh 



20091 ) . Observations of propagating 



disturbances in coronal plumes were interpreted as slow magnetosonic waves not only 
based on their speed but also on their amplitudes' dependence with height, found to be in 
agreement with t he expected variatio n of the wave amplitude with height in gravitationally 



stratified plasma (jOfman et al. 



19991 ). On the other hand, some recent studies suggested 



that propagating EUV disturbances are most likely intermittent transient upflows because 



their association with persistent outflows at the loop base (e.g., 



Hara et al. 



20081 ). and 



their correlation with blueward asymmetries of line profiles using Hinode/EIS observations 



(IMcIntosh &: De Pontieu 



2009 



Tian et al. 



20111 ) . However, these observational s ignatures 



were found to be consistent with slow mode waves as well ( jVerwichte et al. 



20101 ) . 



Our 3D MHD model results of flows and associated waves in hot coronal loops 
show the close connection between the two phenomena and suggest that both may be 
present simultaneously and produced by the same impulsive events at the footpoints (or 
corona/transition region interface) of coronal loops. In particular, we study the effects of 
localized inflow at the lower coronal boundary in a bi-polar magnetic field of hot AR coronal 
loops. The 3D MHD model reproduced the basic physical coupling between steady and 
periodic flows, waves, and loop density and magnetic structures in realistic AR geometry, 
with applications to coronal seismology, extending previous studies that use more simplified 
models and help addressing the controversy of flows vs. waves discussed above. The main 
limitations of our AR model are the use of idealized magnetic structure with dipole field, 
and gravitationally stratified density in the initial state. In addition the isothermal energy 
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equation is used, simplifying the calculations considerably at the expense of calculating the 
effects of heat conduction, heating and cooling of the material injected into the loops, and 
the possible heating of the loops by the energy released in the impulsive events. 

We investigate the effects of steadily injected inflow in short and long magnetic loops 
and the associated oscillations in the initial impulsive stage. We find that the inflow at 
the footpoints of coronal magnetic loops leads to formation of higher density than the 
surroundings loops due to the siphon effect in the curved bipolar magnetic structure. The 
inflows excite coupled longitudinal and transverse oscillations in the loops in the initial 
(impulsive) phase that later damp primarily through leakage as the loops reach the new 
steady state with background flow. The fundamental standing slow mode waves are seen 
to form in short and long loops in agreement with past observations and models discussed 
above. The fast magnetosonic waves induced by the inflow and nonlinearity are detected 
as well. In short loops with steady inflow in low-/? regions of the AR the properties of the 
oscillations (period, wavelength, damping rate), depend on loop length, and curvature of the 
magnetic structur e in qualitative agreement w i th previous studies of waves i n loops with out 



backgr ound flow ([McLaughlin fc Ofman 



2008 



Selwa fc Ofman 



2009 



2010 



Selwa et al. 



2011a 



pj). We also find that the flows produce impulsively nonlinear fast mode waves in the 



loops. 

We investigate the effects of quasi-periodic flow injection in the localized region at 
lower coronal boundary, and find that driven undamped oscillations are produced in the 
magnetic loop of the bipolar AR's. As a result, the average density of the loops increases 
as well in this case. The period, amplitude, and damping rate of the driven oscillations are 
determined primarily by the properties of the driving periodic flows such as their period 
and amplitude, and the nonlinear interaction with the loop structure. In hot and long loops 
the quasi-periodic inflows with subsonic magnitude at the lower boundary can produce 
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shock-like nonlinear compressive disturbances that propagate obliquely to the magnetic 
loop direction. 

We conclude that, a way to determine whether the observed loop oscillations result 
from slow waves or quasi-periodic flows is to evaluate the relation between the oscillation 
properties such as period, phase speed, phase relations, damping time, and the various 
loop parameters such as loop length, density, temperature, and magnetic field structure, 
since these relations are different for waves and for quasi-periodic flows as demonstrated 
in our study. The presence of the temporal damping of the oscillations on the time scale 
shorter than the decay time of the flow (as in the case of steady inflow) supports the 
wave interpretation of the observed oscillations. In the case of periodic inflow studied here 
the oscillations do not exhibit damping. We conclude that the excitation of slow-mode 
(and some transverse) oscillations and flows observed with SUMER, TRACE, Hinode, and 
recently with SDO in coronal loops may result from the impulsive flow injection of plasma 
at the corona- chromosphere interface at the footpoints of the coronal AR loops. Thus, our 
study suggests that flows and waves in coronal loops may result of the same impulsive 
events that drive these phenomena. 
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Fig. 1. — The initial dipole magnetic field (white curves) used for the model AR. The 
intensity scale shows the magnetic field magnitude at the base of the AR and the white 
arrows on some field lines show the direction of the magnetic field. The dimensions are in 
normalized units. The arrow shows the location and direction of the injected flow. 
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Fig. 2. — The initial density (left), fast magnetosonic speed (middle), and plasma /3 (right) 
in the xz plane at t = in the model AR. The contours on Vf show the 50%, 25%, and 
12.5% levels of the maximal value, and the contour (3 = 1 shows the heights where this value 
is reached. 
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Fig. 3. — The initial structure of the long (solid line style) and short (dashes) loops, (a) The 
magnetic field lines at the center of the loops in the xz plane at y — 0. The red circles A, B, 
C, D mark the locations of temporal evolution plots of the variables shown below, (b) The 
magnetic field intensity along the loops, where s is the coordinate along each loop, (c) The 
density along the loops, (d) The fast magnetosonic speed, Vf(s) along the loops. 
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Fig. 4. — The snapshots of the density and the velocity in the xz plane at t — 28OT4 with 
steady inflow (Vo = 0.05). The intensity scale shows the magnitude of the density and the 
velocity, and the arrows indicate the direction of the flow for velocities greater than 6% of the 
maximal value, (a) Density for the short loop; (b) velocity in the short loop; (c) density for 
the long loop; (d) velocity in the long loop. The rectangles marked E-H mark the locations 
of the time-distance plots shown in Figure [5j Animations of the density are available in the 
online version of the journal. 
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Fig. 5. — The time-distance plots for the boxed regions E-G of the short loop and H of the 
long loop shown in Figure HI Panel E is the time-distance plot obtained from the cut E in 
the z direction of the short loop centered at x — 0. Panel F is the time-distance plot of the 
cut F in the x-direction centered at z = 1.82. Panel G is the time-distance plot of the cut G 
in the x-direction centered at z = 1.27 of the short loop. Panel H is the time-distance plot 
of the cut H in the x-direction centered at z — 1.82. In panel H the initial transient time 
is removed to show better the subsequent evolution of the loop's density. In all panels the 
temporal cadence is 5ta- Both, displacement and density oscillations are present showing 
evidence for the transverse and compressional waves in the x — z plane of the loops. 
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Fig. 6. — Same asHJ but in the xy plane for short loop, (a) Density; (b) velocity. Animation 
of the density is available in the online version of the journal. 
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Fig. 7. — The temporal evolution in a short loop of the velocity components (V x : dots; 
V y : dashes; V z \ solid) (top panel), the perturbed magnetic field components (B x : dots; B y : 
dashes; B z : solid; |A_B|: dot-dashes) (middle panel) , and the perturbed density (lower panel) 
at the right footpoint position A = (0.59, —0.01, 1.27) for steady inflow. The red dashed lines 
through the third peak of V z and the peak of \AB\ help evaluate the phase shifts between 
these variables. 
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Fig. 8. — The temporal evolution in a short loop of the velocity components (V x : dots; 
V y \ dashes; V z : solid) (top panel), the perturbed magnetic field components (B x : dots; B y : 
dashes; B z : solid; |A_B|: dot-dashes) (middle panel), and the perturbed density (lower panel) 
at the apex of the loop B = (0.01, —0.01, 1.6) for steady inflow. The red dashed line marks 
three peaks of An helping evaluate the phase shift between the variables. 
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Fig. 9. — The temporal evolution in a long loop of the velocity components (V x : dots; 
V y : dashes; V z : solid) (top panel), the perturbed magnetic field components (B x : dots; B y : 
dashes; B z : solid; \AB\: dot-dashes) (middle panel), and the perturbed density (lower panel) 
at the right footpoint position C = (1.08, —0.01, 1.82) for steady inflow. The red dashed lines 
through the third peak of V z and the peak of |A£>| help evaluate the phase shifts between 
these variables. 
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Fig. 10. — The density (top panels) and the velocity (lower panels) in the xz plane (y = 
—0.0136) at t — 22.5, 105, 188 ta (from left to right) in a short loop injected periodically 
with flow. The arrows show the direction of the flow and the intensity shows the magnitude 
of the flow and the density. Animations of the density are available in the online version of 
the journal. 
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Fig. 11. — The temporal evolution of the velocity components (V x : dots; V y : dashes; V z : 
solid), the perturbed magnetic field components (B x : dots; B y : dashes; B z : solid; |A5|: dot- 
dashes), and the perturbed density at the right footpoint position A = (0.59, —0.01, 1.27) of 
the short loop with periodic inflow. The red dashed through the 4th peak of V z line helps 
evaluate the phase shift between the variables. 
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Fig. 12. — The density (top panels) and the velocity (lower panels) in the xz plane at 
t = 22.5, 105, 188 ta (from left to right) in a long loop injected periodically with flow. The 
arrows show the direction of the flow and the intensity shows the magnitude of the flow and 
the density. Animations of the density are available in the online version of the journal. 
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Fig. 13. — The temporal evolution of the velocity components (V x : dots; V y : dashes; 
V z : solid), the perturbed magnetic field components (B x : dots; B y : dashes; B z : solid; 
\AB\: dot-dashes), and the perturbed density near the apex of the loop at position 
D = (0.01, —0.01, 2.5) of the long loop with periodic inflow. The red dashed line though the 
3rd peak of V x helps evaluate the phase shift between the variables. 



Fig. 14. — Perturbed density isocontours due to inflow in the short loop: (a) steady inflow, 
isocontour level at p = 0.12; (b) periodic inflow, isocontour level at p = 0.1. Same for the 
long loop: (c) steady inflow, isocontour level at p = 0.035; (d) periodic inflow, isocontour 
level at p = 0.05. The animations of the density isocontours for the four cases are available 
in the online version of the journal. 



